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Abstract 

We discuss how the conditions at high redshift differ from those at low 
redshift, and what the impact is on the galaxy population. We focus in 
particular on the role of gaseous dissipation and its impact on sustaining 
high star formation rates as well as on driving star-bursts in mergers. 
Gas accretion onto galaxies at high redshifts occurs on a halo dynamical 
time allowing for very efficiently sustained star formation. In addition 
cold accretion flows are able to drive turbulence in high redshift disks 
at the level observed if at least 20% of the accretion energy is converted 
into random motion in the gaseous disk. In general we find that the 
fraction of gas involved in galaxy mergers is a strong function of time 
and increases with redshift. A model combining the role of dissipation 
during mergers and continued infall of satellite galaxies allows to re- 
produce the observed size- evolution of early-type galaxies with redshift. 
Furthermore we investigate how the evolution of the faint-end of the 
luminosity function can be explained in terms of the evolution of the 
underlying dark matter evolution. 



1 Introduction 

There is ample evidence that the conditions under which galaxies formed at 
high redshift where quite different to those at later times during the evolution 
of the universe. Impressive evidence supporting such a view has been lately 
collected from detailed integral field spectroscopy of z ^ 2 galaxies (Forster 
Schreiber et al. 2006, Genzel et al. 2006). These observations reveal massive, 
high star forming galaxies that show rotational structure that resemble those 
of disk galaxies. In contrast to low redshift disk galaxies however, the ratio of 
gas rotational velocity to velocity dispersion is only of the order few compared 
to V"/cr ~ 10 at low z. Furthermore it has been well estabHshed by now, that 
the star formation rate in the universe is a decHning function of time (Hopkins 
2004) and that this is driven by a down-sizing in the average star formation 
rate of galaxies of a given mass as a function of time (Juneau et al 2005, 
Noeske et al. 2007). Besides star formation, merger rates (e.g. Le Fevre 
et al 2000) and AGN activity (Hasinger et al. 2007) show an upward trend 
with redshift, peaking around z ^ 3 — 4 indicating the importance of the 
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high-z universe for the assembly and formation of galaxies. That indeed a 
large fraction of the stellar mass of massive galaxies is already in place at 
z>2 finds it support from studies of their stellar population (Thomas et al. 
2005) and the evolution of the mass density in units of the z = Q mass density 
(Perez-Gonzalez et al. 2008). 

The hierarchical ACDM paradigm has been successful in predicting and re- 
producing several of the above mentioned trends. In particular the increasing 
merger fraction as a function of redshift is a natural outcome (kb2001, ks2008, 
farouhki ma) of this model and driven by the merging history of the under- 
lying dark matter halos. With respect to the evolution of baryonic physics 
within such evolving haloes early work by (Binney 1977, Silk 1977, Rees & 
Ostriker 1977) laid out the ground work. Radiative cooling and the collapse 
into a rotationahy supported disk of initial hot gas at the halos virial tem- 
perature, followed by star formation in this disk where the main mechanism 
to build up galaxies. Besides its simplifying assumption, e.g. about spherical 
symmetry of the halo, its profile and other properties of the dark matter and 
associated gas infall, these models proved to be very successful in predicting 
galaxy properties and in particular a transition mass scale at which galaxy 
formation becomes very inefficient due to long radiative coohng times of the 
gas (Dekel & Birnboim 2006), thus causing star formation to slow down or 
even halt. In comparison to detailed high resolution numerical simulations 
these models have shown to predict cooling rates too low by a factor of a 
few, especially in low mass halos and at early times (Keres et al 2005). This 
discrepancy can be attributed partly to the geometry of the accretion process 
in cosmological simulations which is not necessarily spherical symmetric as 
assumed in simplified models, but fohows the cosmic filaments. However, the 
essence of the accretion process of cold gas at high redshift can be summarized 
by fast and efficient, on the host halos dynamical time. In effect cold gas is 
accreted with a rate comparable to the dark matter accretion rate times the 
cosmological baryon fraction fbMnM- This situation will not hold in massive 
halos however, where stable shocks heat the infalling gas to the halos virial 
temperature and cooling times are longer than the dynamical time, and nei- 
ther at late times when already a substantial fraction of baryons have made it 
into stars. In the following section we will investigate the impact this efficient 
provision of cold gas has on the evolution and formation of galaxies. 

2 The Model 

We use the semi-analytic model Simurgh to model the formation and evolu- 
tion of galaxies. The dark matter history is calculated using the merger tree 
proposed by Somerville et al. 1999 with a mass resolution of 2 x lO^M©. 
The baryonic physics within these dark matter haloes is calculated following 
recipes presented in Springel et al. (2001) and references therein, including a 
model for the reionizing background by Somerville 2002. In our simulation, we 
assume that ehiptical galaxies form whenever a major merger (M1/M2 < 3.5 
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with Ml > M2) takes place. We assume that during this process, all the cold 
gas in the progenitor disks will be consumed in a central starburst, adding to 
the spheroid mass, and that all stars in the progenitor disks will contribute to 
the spheroid as well. Furthermore, we also add the stars of satellite galaxies 
involved in minor mergers to the spheroid. The merger time scale for galaxies 
is calculated using the dynamical friction prescription in Springel et al. (2001) 
and we find that the predicted merger rate is in good agreement with obser- 
vations Khochfar & Burkert 2001. For more modeling details, we refer the 
reader to Khochfar & Burkert 2005 and Khochfar & Silk (2006a). Please note 
that our simulation does not include AGN-feedback (Schawinski et al. 2006) 
or environmental effects (Khochfar & Ostriker 2007) that have influence on 
the most massive galaxies. Throughout this paper, we use the fohowing set 
of cosmological parameters derived from a combination of the 5-year WMAP 
data with Type la supernovae and measurements of baryon acoustic oscilla- 
tions (Komatsu et al. 2008): f2o = 0.28, Oa = 0.72, Ob/fio = 0.16, erg = 0.8 
and h = 0.7. 

2.1 Evolution of the faint-end luminosity function 

The faint-end of the galaxy luminosity function offers a strong constraint on 
the efficiency of galaxy formation within small mass dark matter halos. Cos- 
mological N-body simulations predict a power-law slope of a ~ —2 at all 
redshift. Observations of the faint-end luminosity function however, show a 
slope that is less steep and varies depending on the observed band and en- 
vironment between a <^ —0.9 to a ^ —1.5 (Ryan et al. 2007). The widely 
accepted explanation is that feedback from supernovae is hindering star for- 
mation in low mass haloes. The rational behind this argument is that the 
specific energy of hot gas is lower in low mass halos and therefor the amount 
of cold gas reheated per supernovae is larger than in massive halos (Dekel & 
Silk 1986). The way that supernovae operate per solar mass stars formed is 
independent of redshift considering the same stellar initial mass function, thus 
it comes as a surprise that observations with the Hubble Space Telescope a 
steepening of the slope with redshift (Ryan et al. 2007). 

In Fig. 1, we show the predicted evolution of a{z) for our best-fit local 
model, i.e. a model that is chosen to best fit the local luminosity function 
(Khochfar et al. 2007). For consistency with the majority of observations, we 
calculate the faint-end slope for the rest-frame FUV at 2: > 4 and at lower 
redshifts for the rest-frame B-band. We indeed find an evolution in a with 
redshift that is in fair agreement with the observed evolution. 

The immediate question that arises is, what infiuences and is the main 
driver for the evolution in al Generally, supernova feedback is considered the 
dominant mechanism in shaping the faint-end of the luminosity function Dekel 
& Silk 1986. The shaded region in Fig. 1 shows the range of linear fits to 
a{z) that we find by varying the star formation efficiency and the supernovae 
feedback efficiency. We infer only a very modest change in a{z) for reasonable 
choices of feedback efficiencies, and therefore conclude that another process 
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Figure 1: Left figure: The slope a at different redshifts as predicted by the 
best fit local model. FiUed symbols show results from the simulation and 
the solid line is the best fit to the simulation data. Errorbars indicate 1 — a 
errors. The dashed lines show the fit to the compiled data in Ryan et al. 
(2007). The shaded region shows the range of linear fits to a{z) that we find 
when varying the star-formation and supernovae feedback efficiency. Stars are 
the compilation from Ryan et al. (2007). Right figure top panel: Relation 
between UV luminosity and host dark matter halo mass for central galaxies 
at ^; = and z = Q. Bottom panel: Dark halo mass function at the same 
redshifts. 



must be responsible for the observed evolution in oi{z). 

The mass function of dark matter halos is known to show a strong evo- 
lution with redshift. The galaxies contributing to the luminosity function 
around L* are mostly central galaxies in their dark matter halos, i.e. the 
most luminous galaxy within the halo (e.g. Khochfar & Ostriker 2008). It 
is therefore not unreasonable to assume a connection between the evolution 
of a{z) and that of the dark matter mass function. When considering the 
luminosity of central galaxies residing in dark matter halos of the same mass 
at different redshifts, we find that at early times, central galaxies are up to 
three magnitudes brighter than their counterparts in low redshift halos (see 
Fig. 1 right upper panel). This is even the case for halos hosting sub-L* 
galaxies. Similar results have been reported by (Kobayashi et al. 2007), who 
showed that dwarf galaxies at early times are not affected by supernova feed- 
back in their simulations because cooHng times are very short in these halos 
and the energy injected by the supernovae is rapidly dissipated away. The 
slope in the region of dark matter halos that host sub-L« galaxies is steeper 
at high redshift, and consequently so is a (Fig. 1 right lower panel). The 
same is true for other choices supernovae feedback and star formation effi- 
ciency, thereby explaining why we do not find any strong dependence of a{z) 
on these parameters. 
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Figure 2: Left figure: Observed and modeled size-evolution of early-type 
galaxies of various masses. Right figure: Size-evolution as a function of stellar 
mass. We find that massive elhptical galaxies show the strongest size-evolution 
with up to a factor of ~ 4 from ^; = 2 to ^; = 0. 

3 Size-evolution of early-type galaxies 

Various observational surveys have revealed that elliptical galaxies are more 
compact at high-z (Trujillo et al. 2006). A natural way to explain such a 
behavior is the combination of dissipation during merger, which drives gas 
to the centre of the remnant and makes it more compact than in the case of 
no dissipation (Cox et al. 2006, Khochfar & Silk 2006b) and the continued 
accretion of satellites that puff up the host galaxy (Naab et al. 2007, Khochfar 
& Silk 2006b). We find that in general the average gas fractions in mergers 
increase from 10% in massive mergers at ^ = to up to 30% at 2; = 4. In 
addition to that we find that massive galaxies have an order of magnitude more 
minor than major mergers. To estimate the size evolution we adopt a model 
in which the relative size of to remnants of the same mass is proportional the 
relative amount of dissipation they encountered during their merging history. 
In this way a remnant of a gas-rich merger is more compact than that of a 
gas poor merger and equally a remnant of many gas-poor satellite mergers is 
less compact than that of a major merger involving gas. 

In Fig 2 we calculate the size evolution for the same redshifts presented 
in Trujillo et al. (2006). The authors took the mean effective radii of the 
ln(i?e) distribution for galaxies above two mass thresholds of 3 x lO^^h^Q 
Mq and 6.6 x lO^^/ifg^ Mq from the SDDS sample of early-type galaxies and 
divided the effective radii of early-type galaxies at higher redshifts by this 
value. After arranging their galaxies in various redshift bins they calculated 
the means of these ratios and presented these values. We here use the same 
method to compare our results to theirs. For both cases of Hmiting masses, 
the agreement is excellent. It appears that the difference in sizes is more 
significant for massive early type galaxies. In the right part of Fig 2 we predict 
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the size-evolution in four different mass ranges based on the relative amount 
of their merger component. While local early-type galaxies between 10^° M© 
and 10^^ Mq are around 1.25 times larger than their counterparts at z = 2, 
local early types with masses larger than 5 x 10^^ Mq will be approximately 
4 times larger than their counterparts at z = 2. This dramatic change in 
sizes in our model results from massive galaxies at high redshifts forming in 
gas-rich mergers while galaxies of the same mass at low redshifts form from 
dry major mergers (Khochfar & Burkert 2003, Naab et al. 2006) and minor 
mergers with small total amounts of cold gas. It is interesting to note that 
models with a characteristic shut-off mass scale for coohng of gas predict dry 
mergers as the main mechanism to grow massive galaxies (Khochfar & Silk 
2008a) and hence imply a large size evolution. 



4 Cold accretion and turbulence 

As mentioned in the introduction, observations reveal a large population of 
massive, high star forming disk galaxies with gas velocity dispersion of the 
order 40 km/s (Genzel et al. 2006). One open question is how these velocity 
dispersion can be driven. Possible explanations that have been put forward 
range from gas-rich galaxies mergers (Robertson & Bullock 2008), tapping 
into gravitational energy of the disk and supernovae feedback (Bournaud & 
Elmegreen 2009, Dekel et al. 2009) to the release of accretion energy by cold 
streams (Khochfar & Silk 2008b). We here investigate latter assumption. 
Numerical simulations by Ocvirk et al. (2008) showed that the fraction of 
cold accreted gas that makes it way down to 0.2 Ryir as a function of dark 
halo mass is roughly constant down to z = 2 and a strong declining function 
of halo mass. We model this behavior by enforcing the simulation results onto 
the gas cooling rates in our model. This way we make sure to be consistent 
with the simulation results. Furthermore, we assume that only cold accreted 
gas is able to deposit a fraction rj of its kinetic energy once it hits the galactic 
disk, while hot gas that radiatively cools down does not contribute. Assuming 
equilibrium between the contributions from accretion and dissipation in the 
disk leads to a simplified expression for the gas velocity dispersion as 



CTg = y 'r]V2j^tdynVvir (1) 

with Mace the cold accretion rate. We here assumed that the energy dissi- 
pation time scale in the galactic disk is proportional to the local dynamical 
time, and that the contribution to turbulence in the inter-stellar-medium from 
supernovae is negligible compared to the contribution from accretion, which 
is the case for high accretion rates that are found at ^ > 2. In the left of Fig 3. 
we show the predicted correlation between the ratio of disk rotational veloc- 
ity and gas velocity V/a and disk rotation V from Khochfar & Silk (2008b). 
Overlaid are various observations at high redshift showing good agreement 
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between a model in which eta = 0.2. This result indicates, that only a smah 
fraction is actually necessary to drive sufficiently high turbulence. The strong 
incline of V/a as a function of F is a direct consequence of the smaller cold 
accretion fraction in massive halos, thus this model suggest that the functional 
form of the cold accretion fraction at high z directly relates to a correlation 
of V/a and V. 




Figure 3: The relation between V/a and V for modeled galaxies (small filled 
circles). The observations are from Genzel et al. (2006) (large filled circle), 
Genzel et al. (2008) (filled stars), Stark et al. (2008) (ffiled triangle) and 
Cresci (2008, ApJ submitted) (filled squares). We here assume that 18% 
of the accretion energy from cold fiows is used to drive turbulence in the 
disk. Note that mergers are excluded as well as satellite galaxies that by 
construction do not accrete cold gas anymore. 



5 Discussion 

In this paper we presented results on the formation of galaxies and their 
properties at high redshift. We investigated the origin of the evolution of the 
faint-end luminosity function, finding that in our model it is mainly driven 
by the evolution of the underlying the dark matter mass function. Low mass 
galaxies at high z occupy a range in dark matter halos that lies on a much 
steeper part of the dark matter mass function than their counterparts at low 
redshifts. This is the main reason for an evolution in the model faint-end slope. 
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Supernovae feedback on the other hand does not introduce an evolving faint- 
end slope, though it will change the slope to become shallower at the roughly 
same rate at all redshifts. An interesting possibility in this respect is the 
choice of IMF. In the case of a top-heavy IMF one would expect more feedback 
contribution. It has been argued, that the discrepancy between observed star 
formation rates and integrated stellar mass density in the universe can be 
reduced by applying a top-heavy IMF at z > 2. Our simulations however, 
show that an increase or decrease by a factor of few will not change the general 
trend of the faint-end slope with redshift, thus not drawing support to the 
notion of an IMF change at high z. 

We continued by focusing on the role that the available fuel for star for- 
mation plays. We showed that the formation of galaxies is more efficient at 
high redshifts, which is a direct result of the fast accretion of cold gas at high 
z. The accretion of cold gas at high z parallels that of the dark matter, and 
hence is much faster than at low redshift, resulting in high star formation 
rates that can be sustained to build up sufficient massive galaxies. 

Further consequences of high gas fractions in galaxies at early times are 
galaxy mergers with large fraction of dissipation. Such mergers have shown to 
be result in very compact remnants. We applied a simplified model that scales 
the sizes of remnants with the amount of available fuel during their merger 
history. Such a model proves to be able to reproduce observed size-evolutions 
of early-type galaxies very well. We find the strongest size-evolution for the 
most massive galaxies > 5 x lO-*^^ M©, making a factor of 4 evolution in 
their sizes from 2; = 2 to ^; = 0. The main cause for this strong trend is 
the occurrence of very gas rich mergers at high z while the local counterparts 
formed rather recent from dry mergers and a large fraction of satellite galaxies 
that help to further puff up the host galaxy. 

The main mode of accretion at high redshifts is by cold flows of material 
coming into the halo along with the dark matter halos. In the case of M* halos 
at high z these cold accretion flows reach the host halo along cosmic fllaments 
and go straight down to the centre of the halo, where they can feed the main 
galaxy on a dynamical time. By the time the cold accretion flows reach the 
main galaxy they will have acquired a substantial amount of kinetic energy 
that will have to be dispersed. Some if it will be radiated away and some 
of it will go into rotational energy of the newly arrived material. Another 
possible option is that some fraction of this energy is used to drive turbulence 
in the gaseous disk. Observations show that indeed high-z disks tend to have 
higher gas velocity dispersions than local disks. Using results from numeri- 
cal simulations on the fraction of cold accreted material as function of halo 
mass one flnds that only 20% of the accretion energy of cold flows is actually 
needed to drive turbulence at the level that it is observed over a wide range 
of galaxy masses and rotational velocities of disks at high-z. Interestingly one 
can recover the steep correlation between V/a and V by just invoking energy 
equihbrium between accretion energy and energy dissipation in the disk, and 
the functional dependence of the fraction of cold accretion as a function of 
halo mass. It will be interesting now to see how the fraction of cold accretion 
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will reflect on Vja for massive galaxies at low z. 

We here presented a number of model results on galaxy properties at high 
redshift, that rehed on the faster assembly of structure, in particular the 
accretion of cold gas onto galaxies. Coming observations will reveal further 
details of the galaxy population at high redshift, that will certainly require 
more detailed modeling and will allow to test current models. 
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